Occupation time distributions for the telegraph process 
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Abstract 

For the one-dimensional telegraph process, we obtain explicit distribution of the occupa- 
tion time of the positive half-line. The long-term limiting distribution is then derived when 
^^ I the initial location of the process is in the range of sub-normal or normal deviations from 

the origin; in the former case, the limit is given by the arcsine law. These limit theorems 
^ . are also extended to the case of more general occupation-type functional. 



a^ 



1. Introduction 



Let B = (Bt, t ^ 0) be a standard Brownian motion on M starting from the origin (Bq = 0), and 
consider the occupation time functional 

-S: ^T-=^ [ H{Bt)dt, r>o, (1.1) 



where H{x) is the Heaviside unit step function (i.e., H{x) = for x ^ and H{x) = 1 for x > 0). 
That is to say, t}T S [0, 1] is the proportion of time spent by the Brownian motion {Bt, ^ i ^ T) 
on the positive half-line. It is well known that the probability distribution of the random variable 
^ I t)T does not depend on T (which is evident from the scaling property of the Brownian motion and 

-^ ' the fact that H{ax) = H{x) for any a > 0) and is given by the classic arcsine law, 

^H . 2 

cn ■ P{f)T ^y] = - arcsin ^, ^ y ^ 1, (1.2) 

f^ ■ with the probability density 

o' 

Pas(y) := ^^=== , 0<y<l. (1.3) 

7rVy(i -y) 

S^ . The beautiful formula (|1.2p dates back about 70 years to P. Levy |Le401 Theoreme 3, pp. 301- 

H \ 302] , who has also proved that the arcsine law (jl.2p is the limit distribution for the relative frequency 

of positive sums among consecutive partial sums of independent symmetric Bernoulli (0-1) random 
variables |Le40[ Corollaire 2, p. 303]. Using the invariance principle, the latter result was extended 
by P. Erdos and M. Kac |EKa47| to the case of sums of arbitrary i.i.d. random variables with zero 
mean and unit variance (cf. |St931 Theorem 4.3.19, p. 236]). More recently, R. Khasminskii |Kh99j 
obtained the limit distribution, as T — t- oo, of more general functionals of the form 

i)T{x;f):=^J f{x + Xt)dt, 

where {Xf, t ^ 0) is a diffusion process on M {Xq = 0) with generator L = — a(x)d^/dx^, and 
/ : R — )• M is a probing function from a suitable class. In particular, the results of |Kh99j imply that 
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if lim2^^_5>-|-oo cl(x) = ao > and / is a bounded piecewise continuous function such that 

Hm - r f{u)du = fi, f+^f-, (1.4) 

x-5>±oo X Jq 

then the distribution of the random variable {\:)j'ix', f) — f^)/{f+ — /_) converges weakly, as T — t- oo, 
to the arcsine law (jl.2p . 

In the present paper, we obtain similar results for the so-called telegraph process defined by 

Xt := Vo I (-1)^" du, t ^ 0, (1.5) 

Jo 

where {Nt, t ^ 0) is a homogeneous Poisson process (with rate A > 0), Vq is a random variable 

with equiprobable values ±c independent of the process Nf, and c > is a parameter (see |Go5H 

IKa741 IPi91j ). That is, Xt is the position at time t ^ of a particle starting at t = from the 

origin and moving on the line with alternating velocities itc, reversing the direction of motion at 

each jump instant of the Poisson process Nt; the initial (random) direction is decided by the sign 

of Vq. Note that the process Xt itself is non-Markovian, however if Vt = dXt/dt = (— l)^Vo is 

the corresponding velocity process, then the joint process {Xt,Vt) is Markov on the state space 

M X {— c, +c} (see |EtK861 §12.1, p. 469]). We shall also consider the conditional telegraph processes 

obtained from Xt by conditioning on Vq, 



Xf := ±c / (-1)^" du, t ^ 0, (1.6) 

Jo 

where the choice of the + or — sign determines the initial direction of motion. 

Remark 1.1. Here and throughout the paper, we adopt a notational convention that any formula 
involving the zb and =F signs combines the two cases corresponding to the choice of either the upper 
or lower sign, respectively. 

Remark 1.2. The telegraph process is the simplest example of so-called random, evolutions (see, e.g., 
|EtK86[ Ch. 12] and [Pl9Tl Ch. 2]). 

The model of non-interacting particles moving in one dimension with alternating velocities (up- 
dated at random on a discrete time grid) was first introduced in 1922 by G. I. Taylor |Ta22j 
in an attempt to describe turbulent diffusion; later on (around 1938-1939) it was studied at 
length by S. Goldstein |Go51j in connection with a certain hyperbolic partial differential equa- 
tion (called the telegraph, or damped wave equation, see (|1.7p below) describing the spatio-temporal 
dynamics of the potential in a transmitting cable (without leakage) |We55j . In his 1956 lecture 
notes, M. Kac (see |Ka74] ) considered a continuous-time version of the telegraph model. Since 
then, the telegraph process and its many generalizations have been studied in great detail (see, 
e.g., |Or90t IPi91t IOr95t IRa991 IWe02j ). with numerous applications in physics |We02] . biology 
|Ha991 IHH05J . ecology |OL01j and, more recently, financial market modelling [Ra071 IRM08J (see 
also further bibliography in these papers). 

An efficient conventional approach to the analytical study of the telegraph process, analogous 
to that for diffusion processes, is based on pursuing a fundamental link relating various expected 
values of the process with initial value and/or boundary value problems for certain partial differential 
equations (see, e.g., |Go51t [Or901 IOr9"5l IRa971 IRa99( IRa06j ). In particular, Kac |Ka74] has shown 
that, for any bounded continuously differentiable function (70 ^ I^ ~^ I^i the functions 

v^{x,t) :=E[go{x + Xf)], x e R, t ^ 0, 

satisfy the set of partial differential equations 

dv^(x,t) dv^ix,t) , / _i_, X , x\ 

k^^"" ^y ^ =TA(^+(x,t)-7;-(rr,t)), t > 0, 
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with the initial conditions 

v^{x,0) = go{x), X €R. 

These equations can be easily combined (see details in [Ka74| or |EtK86( §12.1, p. 470]) to show 
that the function 

v{x,t) ■.= E[go{x + Xt)] = -v~{x,t) + -v+{x,t) 

satisfies the telegraph (or telegrapher's) equation (see, e.g., |We551 §15]) 

d'^v , dv n d'^v , „, 

with the initial conditions 

dv 
^(x,0)=5o(x), — (x,0) = 0. (1.8) 

Remark 1.3. The telegraph equation (jl.Zp first appeared more than 150 years ago in work by 
W. Thomson (Lord Kelvin) on the transatlantic cable [Th54| . 

The (unique) solution of the Cauchy problem ()1.7p - (|1.8p can be written explicitly (see, e.g., 
|We55[ §§46, 74] or (PMI §0.4]) as 

v{x,t) = -e-^\gQ{x + ct)+gQ{x-ct)) 



2 

rt 



+ ^e~^' J ^go{x + cu)(xio{XVt^-u^) + —^^h{XVt^-u^))du, 



iU 



where 



n=0 ^ ^ n=0 ^ ' 

are the modified Bessel functions of the first kind (of orders and 1), respectively) |AS721 9.6.12, 
p. 375; 9.6.27, p. 376]. 

It is well known that, under a suitable scaling, the telegraph process satisfies a functional central 
limit theorem. 

Theorem 1.1. Assume that A, c — t- +oo in such a way that c^/A — t- 1. Then the distribution of 
the telegraph processes [X^ , i ^ 0) converges weakly in C[0,oo) to the distribution of a standard 
Brownian motion (Bt, t ^ 0). The same is true for the unconditional telegraph process {Xt, t ^ 0). 

As was observed by Kac |Ka741 p. 501], this result formally follows from the telegraph equation 
(|1.7p , which in the limit A, c — )• +oo, c^/A — )• 1 yields the diffusion (heat) equation 

dv 1 d'^v 

at ~ 2 a^' 

associated with the standard Brownian motion Bf. A rigorous proof of Theorem 11.11 along with 
some extensions, can be found in |EtK86l §12.1, p. 471] and |Ra991 Theorem 5.1]. 

Our main goal in the present paper is to analyze the distribution of the occupation time of 
the telegraph process Xt and, in particular, to obtain a limit distribution, as T — )• oo, of the 
occupation-type functionals of the form riT{x; f) := T^^ L f{x + Xt) dt for a suitable class of 
probing functions /. In particular, we prove that the limit distribution is given by Levy's arcsine 
law providing that the starting point x is in the range of subnormal deviation from the origin (i.e., 
X = o{\/T)). For technical simplicity, we impose a stronger condition on the asymptotics of / at 
iboo, assuming that the corresponding limits f± exist. 
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The rest of the paper is organized as fohows. In Section [2] we state the main results of this 
work (Theorems I2.1h - |2^ . which are then proved in Sections HH?! respectively. Section [3] contains 
a suitable version of the Feynman-Kac formula, with applications to the Laplace transforms for 
the occupation-type functionals under study, which is instrumental for our techniques. We finish 
in Section [8] with concluding remarks and some conjectures, which are illustrated by the results of 
computer simulations. Appendices lAl and IBJ contain alternative (probabilistic) proofs of Theorems 
12.21 and 12.31 respectively. 

2. Statement of the main results 

For T > 0, x € M, consider the following occupation time random variables 

^t{x):=^J Hix + Xt)dt, 4{x):=^j H{x + Xt)dt, (2.1) 

where H{x) = l(o,oo)(^) is the Heaviside step function and Xt, X^ are the telegraph processes 
introduced above (see (jl.Sp . (jl.6p ). Note that the total time spent by the processes {x + X^ , ^ 
t ^ T) at the origin almost surely (a.s.) equals zero, since by Fubini's theorem we have 

E l{oy{x + Xf) dt = F{X^ = -x} dt = . (2.2) 

Jo Jo 

Hence, the complementary quantity 1 — f?y(x) a.s. equals the proportion of time spent by the 
processes {x + X^ , ^ t ^ T) on the negative side of the axis, 

1 f^ 
l_7y±(x) = -/ l(-^fi){x + Xf) dt (a.s.), 

and by symmetry (with respect to simultaneous transformations x i— ?• — x, it i— t- =f) it follows that 

'i]^{x) = l-'r]^{-x), xeR. (2.3) 

Let us consider the function iprit) {t ^ 0) defined by 

After the substitution u = ty, we have in the limit as t J, 0, 

1 /■! 1 _ Q-2\Tty 1 /-i 1 1 

(see (|1.3p ). and so (pxi') is continuous at zero (and hence everywhere on [0, oo)). Note the following 
useful scaling relation, which easily follows from the representation of if given by (|2.5p : 

ipo^rit) = ^T{at), t^O, a > 0. (2.6) 

Let us also set 

i;T{y)--=2XT^T{y)M^-y), O^y^l. (2.7) 

We are now ready to state our first result. 

Theorem 2.1. The random variables rjrp^O) defined in (|2.ip have the distribution 

F{4{0)£dy}=2ipTil)5,±{dy) + ijTiy)dy, sj y ^ 1, (2.8) 

where 5x is the Dirac measure {of unit mass) at point x, with x" = 0, x"*" = 1. Furthermore, the 
distribution of r]T{0) (see (j2.ip ) is given by the formula 

¥{rjT{0)edy}=^T{l)6o{dy) + ^T{l)Si{dy) + ijT{y)dy, ^ y ^ 1. (2.9) 

4 
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In other words, the distribution of i]^{0), r]j^{0) has a discrete part with atom of mass 2(/?t(1) 
at point or 1, respectively, and an absolutely continuous part with the density tpT defined by 
(|2.7p . Similarly, the distribution of ryr(O) has atoms at points and 1, both of mass (/'^(l), and an 
absolutely continuous part with the density ^l^T as above. 

Remark 2.1. The it-duality in ()2.8p becomes clear from relation ()2.3p (with x = 0) and the symmetry 
property il^riy) = '4't{'^ -y) (see ^1^ ). 

Remark 2.2. Using an integral formula (see |AS72l 9.6.16, p. 376]) for the modified Bessel function 
Iq, it is easy to check that the function ^pT defined by (j2.4p admits another representation, 

rXTt 



^T{t) = i^yflj e-yio{y)dy, t > 0, 



which is further evaluated (see [AS721 11.3.12, p. 483]) to yield iprit) = ^e^^^*(/o(ATi) + /i(Art)). 
Thus, the distribution of ??^(0) and r/r(0) can be expressed through the modified Bessel functions 
Iq and Ii, as well as the distribution of the telegraph process (cf. ()1.9p l. 

In the next theorem, we give explicit integral formulas for the distribution of rij,{x) in the case 
X 7^ 0. For simplicity, we only present the answer for x < 0, the case a; > readily following in view 
of the duality relation (j2.3p . 



Theorem 2.2. Assume that x < and set Tq := \x\/c. Then, for any T > 0, the random variables 
rjrp{x) defined in (j2.ip iiave the following distribution: 

(a) if T ^ To then F{r]^{x) = 0} = 1; 

(b) if T>Tq then, for < y ^ 1 - Tq/T, 

¥{r,${x) e dy} = (J^ QtAu)du\ 6o{dy) + 4{dy) + ^^iy,T)dy, (2.10) 

where fi^{dy) := and 

4{dy) := 2e-^^Vr(l - n/T) S^^To/Hdy) + e'^^" ^^t-t, {y^JTJ^) T^kjT ' ^^'^^^ 

r{^-y)T / V \ du 

^t{y, T) := 2rQ±,((i - y)T) My) + J^ Q-Au) ^T-u [j^;;^) y^^ ' (2-i2) 

with ifT and iI^t given by ()2.4p and (j2.7p . respectively, and the functions Q__^{u) (— x > 0) defined 
for all u £ [Tq, oo) by 

\rp p— Am 

QUn):= ,f /i(A7^?^7f), (2.13) 

V" - -'o 

QZM ■■= Ae-^"/o(Ax/^?^lf ) - ^^^^^^l^/i(Ay^?^Tf ). (2.14) 

For the next theorem, we need a few notations. For a > 0, consider the function 

a fa? 



with Laplace transform (see |AS72l 29.3.82, p. 1026]) 

1 

e-'^qa{t)dt = e-''^' , s ^ 0. (2.16) 

5 



oo 



where 



m, 
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Let Ya (a ^ 0) be a family of random variables with values in [0, 1], such that Yq has the arcsine 
distribution ()1.2p . with the density pas (see ()1.3p ). while for a > the distribution of Ya is given by 

nYa e dy} = ma6o{dy) + fa{y)dy, (2.17) 

a := / qa{u) du = ^= e~y /^ dy, (2.18) 

Ji \/2tt Jo 

fa{y)-= -. Pasi-. dn=^===/ ^du. 2.19 

Remark 2.3. It is easy to verify, either from ()2.15p or using the Laplace transform ()2.16p . that 

9a — ?• (^o as a — 7- 0+, where 5o{-) is the Dirac delta function and — t- denotes weak-* convergence of 
generalized functions; hence nia — ?• (which can also be seen directly from the right-hand side of 

(|2.18p ) and fa — )■ Pas (see the first part of formula (|2.19p ). That is, Ya ^ Yq as a ^ 0+, and so the 
distribution of Ya is continuous in parameter a £ [0, oo). 

Theorem 2.3. Suppose that the initial position X^ = x, as weU as the parameters c and A, may 
depend on T in such a way that XT — t- oo and {(P'T/X)~^''^ x— T-aGMasT— ?• oo. Then, as T — t- oo, 

ryHx), r?^(x) A ( f- "5°' (2.20) 

y l-Ya, a^O. 

In particular, for a = the limit is given hy the arcsine distribution (|1.2p . 

To order to generalize these results in the spirit of |Kh99] . let / : M — ^ M be a bounded, piecewise 
continuous function (i.e., continuous on M outside a finite set Df, where it has finite left and right 
limits), such that, for some finite constants f^ ^ f- , 



hm /(x) = /_, hrn /(x) = /+ . (2.21) 

X— >— oo X— >+oo 



Consider the random variables 



^^Ti^-J) ■■= JJ^^J^ I {fix + Xt)-f-)dt, xG 



(2.22) 



Clearly, by a linear transformation of the function, 

fix) - /_ 



fix) ^ fix) :-- 



f+-f^ 

we may and will assume without loss of generality that /_ 

rT 



4i^-J)--=^l fix + Xt)dt. 



X G R, 










(2.23) 


0, /+ = 


1, 


so 


that 


(|2.22p is reduced to 


~-)dt. 










(2.24) 



Theorem 2.4. Let the function f satisfy the above conditions including assumption (j2.2ip with 
/_ = 0, /+ = 1. Suppose that the hypotheses of Theorem 12.31 are satisfied, and assume in addition 
that c^T/X -^ oo as T ^ oo. Then the distribution of r]rpix; f) converges weakly, as T — t- oo, to 
the law determined by the right-hand side of ()2.20p . 

Remark 2.4. Theorem 12.31 may be inferred from the diffusion approximation (Theorem I l.ip of the 
telegraph process (see an alternative proof in the Appendix|B]). However, we will give a direct proof 
of Theorem 12. 3^ which may be of interest in its own right and will also be instrumental for laying 
out the necessary techniques for the proof of Theorem 12.41 where the "diffusion approximation" 
trick does not seem to be readily applicable. 

6 
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3. The Feynman Kac formula and applications 

Let us recall the Feynman-Kac formula for the telegraph processes. 

Theorem 3.1. Let {X^ , t ^ 0) be the telegraph processes (|1.6p . Suppose that qq and g are bounded 
functions on M such that go G C^(M) and g is piecewise continuous, i.e., g S C{M.\Dg), where Dg is 
a finite set, and moreover, f has finite left and right limits at the points of Dg . Then the functions 



v'^{x,t) :=E 



5o(x + X±)exp| / g{x + X^)du 



x£R, t ^ 0, (3.1) 



for all {x, t) gR X IR+ such that x ziz at ^ Dg satisfy the set of partial differential equations 

dv^(x,t) dv^(x,t) , / _i_/ N / x\ / N +/ N 
it ^ "" 'k^ " ^ K(^>0 - v~{x,t)) +g{x)v^{x,t), 

with the initial conditions 

v^{x,0) = go{x), X e M. 

This theorem is proved (see details in |Ra06j ) similarly to the analogous result for diffusion 
processes (cf., e.g., |IM741 §2.6]). An alternative probabilistic representation for the solution of a 
deterministic telegraph-like equation is developed in JDMT08] . 

Let r]^{x) be defined by dHJ. For /3 G M, set 

4(C,t) := E[e-^*''Tt(^rO]^ ^^-^^ i ^ 0, (3.2) 

or more explicitly (cf. (j3.1|) ) 



4(e,t)=E 



expj^^ Hicn + Xt)du^ 



C G M, t^O. (3.3) 



Since H{-) is a bounded function, the expectation in (13. 3p is finite for all /3 G M. 
Let us record some simple properties of the function Vj, . 

Lemma 3.2. For each /3 G M and any T > 0, the functions Vrp{^,t) are continuous on R x M_|_ and 

lim v${C,t) = l, lim t;±(^, t) = e"'^*. (3.4) 

Proof. Continuity in t G M+ is obvious. As mentioned above (see ()2.2p ). for any ^o G IK we have a.s. 
that cT^Q + X^ 7^ for all u G [0, Tt] except on a (random) set of Lebesgue measure zero. Since the 
function H is continuous outside zero, this implies that, for such u, H{cT^ + X^) —^ H{cT^q + X^) 
as ,^ — ;• ^0 and hence, by Lebesgue's dominated convergence theorem, f^ H{cT£^ + X^)du — ^ 
Jq H{cT^o + X^)du as ^ — )■ ^o- The continuity of Vrp{-,t) at point .^o now follows by Lebesgue's 
dominated convergence theorem applied to the expectation (j3.3p . since everything is bounded (for 
a fixed t). 

To prove (|3.4p . note that, for T > and each u ^ 0, we have cT$^ + X^ —^ zboo as ^ — )■ ±00. 
Since H is bounded on M, the claim now follows by dominated convergence. D 

From the definition (13. 3p . it is clear that if /3 ^ then, for each ^ G M, the functions Vj^{^, •) are 
bounded on [0,oo), so we can define the Laplace transform 

I" 00 
w^{^,s):= e-'^v^{S„t)dt (s > 0). (3.5) 

Jo 

Lemma 3.3. Set s := s + (3. For any fixed s > 0, the functions Wj. = Wrp{^, s) defined by p.Sp are 
continuous in ,^ G M and satisfy the following set of differential equations 

-^ = XT{w+ - w;^) ±{s + H{cTO) w^^l, C 7^ 0. (3.6) 
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Moreover, 

lim Wj.{—£,,s) = s~ , lim Wj,{^,s) = s~ . (3.7) 

Proof. The continuity of the functions w^{^,s) in ^ follows from the definition ()3.5|) and the first 
part of Lemma |3.2[ Further, applying Theorem 13.11 (with go{x) = 1 and g{x) = —PT~^H{x)), we 
see that the functions v^ = Vrp{^,t) defined by (|3.2|) satisfy the initial value problem 

T-^ + -Q^ = XT{v+-v^)±l3H{cn)v^, t>0, e±i/0, (3.8) 

v${C,0) = l, CgM. (3.9) 

Integrating by parts and using the initial condition (j3.9p . we have 

r e-'' ^""^5'^^ d* = -4 (^' 0) + ^ /°° e"** 4(^, t) dt = -1 + st/;^(C, s). (3.10) 

Applying the Laplace transform (with respect to t) to equation ()3.8p and taking into account (j3.10p . 
we immediately obtain the differential equation (j3.6p . Finally, the boundary conditions (|3.7p readily 
follow from (j3.4p by Lebesgue's dominated convergence theorem applied to (j3.5p . D 



Let us also make similar preparations for the random variables rjrp^x;/) defined in (|2.22p . As 
explained in Section [2] (see (j2.23p ). without loss of generality this definition can be simplified to the 
form (I223D. Consider the functions (cf. dSJ])) 

4(e,t;/):=E[exp(-/3t7?±(cre;/))], ^eR, t ^ 0, 

and the corresponding Laplace transform 



POO 

w^{ts;f):= e-''v${C,t-f)dt, s>0. 

Jo 



Then, again applying Theorem 13.11 (with go{x) = 1 and g{x) = —j3T~^f{x)), similarly to Lemmas 
3.21 and [3.31 one can show that t(;,p = Wrp{^,s; f), for each s > 0, is a continuous bounded function 
of ^ G M, satisfying the differential equation (cf. ()3.6p ) 



"'^ - Ar(w;+ - w") ± (s + /3/(cT^)) t/;^ =F 1, £,£R\Df, (3.11) 



with the same boundary conditions at ±00 as ([37 



lim Wrp{^,s; f) = s , lim Wj,{£^,s; f) = s . (3.12) 

4. Proof of Theorem 12.11 

In what follows, the prime ' denotes the transposition of vectors. Introducing the vector notations 

WT{C,s):={w+iC,s),w^{ts)y, 1:=(1,1)', i:=(l,-l)', 

we can write down equations (|3.6p and (|3.7p in the matrix form. 



d^ 



At{^,s)wt{C,s)-1 (e/0), (4.1) 



lim wt{C,s) = s~^l, lim wt{C,s) = s~^l, (4.2) 

§— >— 00 §— >+oo 
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where s = s + /3 (see Lemma 13. 3|) and 

f XTJi + sJ2 =: At = At{s), C < 0, 
^rii..):=>-TJ..(s^mmV. = [^^,^^^^,,^^^^^^^^^\l ^^„_ (4.3) 

Note that 

Jii = o, Jii = 2-i, J2i = i, J2i = i, (4.5) 

where := (0,0)'. Hence (see (|0]) ) 

AT(s)l = si, AT(s)i = (s + 2Ar)l. (4.6) 

Let us set 

k = k{s) := ^/s{s + 2XT), k := k{s) = y^s{s + 2\T). (4.7) 

Using formulas ()4.6p and ()4.7p . it is easy to check that the matrix At{s) has the eigenvalues ibK(s) 
with the corresponding eigenvectors 

a± = a±{s) := iLkI + si, ATa± = ±Ka± . (4.8) 

In particular, relations (|4.8|) imply that the exponential of At{s) can be represented as follows: 

e^Ti = \ e< (/ + ^c^AT) + \ e'^^ (/ - k-^At) , (4.9) 

where / is the identity matrix. 

Recall that we are looking for a solution to the boundary value problem (|4.ip - ()4.2p continuous 
at the origin. The following lemma gives an explicit form of such a solution. 



Lemma 4.1. For each s > 0, the differential equation (|4.ip subject to the boundary conditions (|4.2 
has the unique continuous solution given by 



WTiS,,s) 

In particular, 



sk + SK 

e-^^— (Kl-si) + s-4, ^^0. 

SK + SK 



(4.10) 



wr{0,s)=^-^±^^^^^^^^. (4.11) 

sk + SK 



Proof. Observe that the step function w'!^{^,s) := (s + /3H{cT^)^ 1 is a particular solution of 
the equation (j4.ip for each s > and all ^ 7^ 0. Indeed, the function w^(-, s) is piecewise constant 
outside zero, hence {d/d^)w'^{^,s) = (^ / 0), whereas, due to ()4.3p and ()4.5p . 

AT{ts)w*Ti^,s) = \T{s + PH{cT0)'^Jil + J2l = l- 

Therefore, a general solution of the linear differential equation ()4.ip can be represented in the 
form (see (fO]) ) 

f e^^«c(s) + s-^l, C<0, 
n,T{ts) = l , ' , (4.12) 

with arbitrary vectors c{s), c{s) (which may also depend on T). Let us now find suitable c{s) and 
c{s) so that the solution wt{-,s) would satisfy the required boundary conditions at infinity and 
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the continuity condition at zero. From the representation (j4.12p it is clear that conditions (j4.2p are 
satisfied if and only if 

lim e^'r^c(s) = 0, lim e^^^c(s) = 0. (4.13) 

§— >— oo 5— >+oo 

Recalling that At{s) = At{s) and using the exponential formula (|4.9p . it is easy to see that 
conditions (j4.13p are reduced to the equations 

(/ - k-Ut) c{s) = , (/ + k-'^Ar) c{s) = , 

which implies that c(s) and c(s) are eigenvectors of the matrices At and At, respectively, with 
the corresponding eigenvalues k and —k. On account of formulas (j4.8p . this immediately gives 
c{s) = C(s)a_|_, c(s) = C{s)d^, with some real- valued functions C{s), C{s). Therefore, after the 
substitution of expressions (|4.8p . formula (j4.12p takes the form 

WT{ts) = < -,~ ~ ' (4.14) 

Furthermore, taking into account the continuity of wt{-,s) at zero, from ()4.14p we have 

C{s) (k1 + si) + s-^1 = C{s) {-kl + si) + s-^1, 
whence, by equating the coefficients of 1 and i on the left- and right-hand sides, we obtain 

J C{s) K + s~^ = -C{s) k + s~\ 

\c{s)s = C{s)s. 
Solving this system of equations we find 

C(s) = Z^""'! , C{s)- ~^'~^ 



sk + SK sk + SK 

and the substitution of these expression into ()4.14p yields the required formula (j4.10p . 

Finally, the expression (14. lip for wt{Q,s) follows from ()4.10p by setting ^ = and using that 
(3 = s — s (see Lemma 13. Sp . This completes the proof of Lemma 14.11 D 

Lemma 4.2. The components Wrp{0,s) (see (|4.1ip ) are explicitly given by the expressions 



2XT 

k+ s {k + s){k + s) 



wj;(0,s) = ^^^+ ^ , ,^^_ , ^, , (4.15) 



2 2XT 

^^(0,s) = ^ + -^-— — -. (4.16) 

K + S {K + S){k + S) 

Proof. From the vector expression ()4.1ip we have 

Note that, according to (|4.7p . 

K^ -s^ = 2XTs, k^ -S"^ = 2XTs, (4.18) 

hence the expression (j4.17p may be rewritten as 

w^{0,s) = — — + 



Sk + SK \ K^ s k± s 

2XT ( s s \ 2XT ,^ , 



SK -|- SK \ «; =F s k±sj {Hi^s){k±s) 
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which is equivalent to (j4.15p . (j4.16p : for instance, for Wq^{0,s) (corresponding to the choice of the 
upper sign in it and =p), from formula (j4.19p we obtain, again using (j4.18p . 

2XT 2XT{k + s) k + s 



wj^{0, s) 



(k — s)(k + s) (k^ — s^)(fi; + s) s{k + s) 
2 K-s 2 2XT 



k + s s{k + s) K + s (k; + s)(k + s)' 
in agreement with (I4.15p . Thus, Lemma 14.21 is proved. D 

Lemma 4.3. Let the function ipT{t) be defined by (j2.4p . Then, for each s > 0, 



/"OO 1 /'OO 1 

/ e-"VT(t)dt = , / e-^*e-''VT(t)dt = 

Jo /« + s ' Jo ^' t^ + 



(4.20) 



and 

[.-' [j\-f'y^r{y)^it - y)iy) it = (^^^;(,^,) . (4.21) 

where s = s + P and k = k,{s), R = k{s) are dehned in (j4.7p . 



Proof. Inserting (j2.4p and changing the order of integration, we obtain 

1 foo -I _ —2XTu / foo „~st 



r(i) 



/■oo 

/ e-^" (1 - e-2^^«) n-3/2 dn. (4.22) 

Jo 



AttXT^s Jo 
Integration by parts via m~^' ^ du = —2 d{u~^''^) yields the right-hand side of (j4.22p in the form 

2XTy/¥s Jo ^ ' 2ArV^ K + s 

and the first formula in (I4.20p is proved. The second one readily follows by the shift s = s + (3. 

Furthermore, using the convolution property of the Laplace transform, the left-hand side of 
()4.2ip is reduced to the product 

/"oo poo -I 

/ e-^*e-/^VT(t)dt X / e-^VT(t)dt = -———-, 
Jo Jo {k + s){k + s) 

according to ()4.20p . which completes the proof of the lemma. D 

Combining Lemmas 14.21 and 14.31 and using the uniqueness theorem for the Laplace transform 
(13.51). we obtain 



4(0, i) = (1 + e^^* + 1 ± e-^*) ^T{t) + 2XT [ e~^y^T{y)Mt " v) dy- 

Jo 

In particular, setting t = \ (see (j3.2p ) and recalling the definition (j2.7p of the function ij^t, we get 

E[e-/3'?T(o)] = (1 + e"/^ + 1 ± e"^) ^t(I) + / e'^^^riv) dy, 

Jo 

and it is evident (in view of the uniqueness theorem for Laplace transform) that the distribution of 

•qrp (0) is given by formula (|2.8p . 

Finally, the result (j2.9p for ?7r(0) readily follows from (|2.8p and the decomposition 

P{r/r(0)Gdy} = ^P{7?+(0)Gdy} + ^P{r/^(0)Gdy} (0 < y «: 1). (4.23) 

Thus, the proof of Theorem 12. II is completed. 
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5. Proof of Theorem [272] 

The plan of the proof below is to calculate the Laplace transform (see (I3.2p and (13. 5p ) 

l-CO 

uj${C,s)= e"^*E[e-^*''T*(^^?)]dt (5.1) 

Jo 

from the explicit (hypothetical) distribution of r]rp{x) given by formula ()2.10p . and to verify that 
the result coincides with formulas (|4.10p obtained in Lemma 14.101 The claim of Theorem 12.21 will 
then follow by the uniqueness theorem for Laplace transform. To be specific, we will focus on the 
w^ case, the proof for w^ being similar. 

Remark 5.1. In Appendix |B] we will give an alternative proof based on probabilistic arguments, by 
reducing the general case rjj. (x) to r]j, (0) via conditioning on the hitting time of the origin. That 
proof will explain how formulas (|2.1Up can be derived (rather than verified); however, in so doing 
the prior knowledge of the distribution of rjrp (0) (provided by Theorem 12. ip is essential. 

Due to the space-time change {x,T) i— )• {cT^,Tt) used in ()5.ip . the time threshold To = \x\/c 
becomes Tq = T|^|, whereas the former condition T > Tq is converted into t > |^|. As a first step 
in the proof, using the probability distribution proposed by the theorem (including its part (a)) we 
can represent the Laplace transform of tr]^^{cT^) as 

where ^Tj. (^, t) (i = 1, . . . , 5) arise from the three parts on the right-hand side of the representation 
(|2.10p . with the last two further subdivided each into two terms, according to (|2.1ip and (j2.12p . 
More precisely, using the scaling property (j2.6p of the function (px and making the substitutions 
y h^ ty and u i— )■ T{u + |^|) wherever appropriate, the functions Jj^ {^,t) can be expressed as 

/■oo 

J^'\C,t) := T J g+,|^|(Tn)dn, (5.3) 

j¥Htt + \m:=2e-^'-'^\^\Mt), (5.4) 

4'^(e,t+|ei):=e"^^l^l^'e"^^^T*(f)^, (5.5) 



"^ 

4'^ (e, t + \C\):=T 1^ e~^y [f'^Qtm {Tin + ICD) ^T(*-«) (^) ^) dy. (5.7) 

Consequently, from (15. 1|) and (15. 2|) we get 

Tz;+(e,s) = -(l-e-l«l^)+^/ e"^*J-«(e,t)dt. (5.8) 



rii), 



Let us now calculate the Laplace transform (with respect to t) of each of the terms J^j, (^, t) 
(i = 1, . . . , 5). In so doing, the next formula will be useful, 

T e-^*Q+^|^|(rt)dt = e«'^-e(^+^^)« (C < 0), (5.9) 
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where At = y^s{s + 2AT) (see (j4.7p ). which immediately follows from the definition ()2.13p according 
to [XS721 29.3.96, p. 1027]. 

Remark 5.2. An analogous formula for Q~ (needed for the proof in the case of w~) follows from 
(l2l4D by applying [M72l 29.3.93, 29.3.96, p. 1027]). 



(i) From (15. 3p we obtain, integrating by parts and using (|5.9p . 

/ e-''ji}\tt)dt = Ts-^e~'\i\ Q+y|^|(Tt) dt - Ts^W e"^*Q+.|^|(rt) dt 

= 1 (e^« - e(«+^^)«) - - (e^^ - e(''+^^)«) = i (e'^ - e<). (5.10) 

s s s 



(ii) After the substitution 1 1— )• t + |^|, from (|5.4p we get, using formula (|4.20p in Lemma [ 

/ e-^*4')(C,t)di = 2e(^+^^)« / e-('*+«VT(t)dt = 2e(^+^^)«^. (5.11) 

J\i\ Jo ^ + ■5 

(iii) Likewise, from ()5.5p we obtain, recalling the definition ()2.7p of the function ipT and again 
using the scaling property (|2.6p . 

'1^1 Jo \Jo \tJ t J 

roo / rt 



poo /"oo / nt 

/ e-^*4')(e,t)dt = e(^+^^)« / e-^M/ 
ilfl Jo \Jo 

= e{s+\T)i 2XT r e^^* ( f e-^^ friy) fT{t - y) dy") dt 



{k + s){k + s)' ^ ' 

as follows from formula ()4.2ip in Lemma |4.3[ 

(iv) Similarly, taking advantage of the convolution theorem, the Laplace transform of ()5.6p can 
be written as 

I" e-^* 4") (^, t) dt = j^ e-^* (2T 1^ e-^y Mv) Qtm (^(* + 1^1 - y)) dy) di 

/"OO /"OO 

= 2y e-^*e-^VT(t)di xe^^ry e-^*Q+^|^|((t + |e|)r) dt 

), (5.13) 



2 



fe^ - e('*+^^)5 



K + S 

according to formulas ()4.20p and (j5.9p . 

(v) Interchanging the integrations, we can rewrite (|5.7p in the form 

4'^(e,i + 1^1) = Ty|*Q:^|^|(r(n+ lel)) (^*~"e-^^^T(t-n) (^) ^) dn, 

hence, by the convolution theorem, the Laplace transform of J'j^ i^i'^) is reduced to 

y" e-^* 4'^ (e, t) dt = e^« T j^ e-^* Q+^|^| {T{t +\i\)) dt x j^ e"^* (y|* e"^^ V'Tt (f ) y ) di 

= (e< - ^is+XT)^) '^-^^— , (5.14) 

as was shown in (|5.12p and (j5.13p . 

13 



Leonid Bogachev and Nikita Ratanov 

Finally, substituting the results (lETOD . (|5lH) . (l5J2]l . (l5J3]l and (l5H]l into formula ([53]) and 
recalling the expressions (|4.15p . (|4.16p for Wj,{0,s), we get 

Kic, s) = -{i- e<) + e< f ^ + rrrw^-^ 

= - (1 - e'^^) + e'^^ u;+(0, s) = &< fu;+(0, s) - -^ + - , 
s ' \ s J s 

which is consistent with the expression ()4.10p for w'^{^,s) obtained in Lemma |4.1[ Thus, the proof 
of Theorem 12.21 is completed. 

6. Proof of Theorem [231 

It suffices to prove the theorem for the conditional versions ry^, {x) only; indeed, since the latter have 
the same distributional limit, the result for rjT{x) will readily follow (cf. (j4.23p ). 

In the next lemma, we find the Laplace transform for a suitable parametric family Ya (t) extending 
the random variables Ya introduced in Section [2] (see (j2.17p . (|2.18p and (j2.19p l. Recall that s = s+/3. 



Lemma 6.1. For any a ^ and t > 0, set Ya{t) := tY^in. Then, for any s > and /? > 0, we have 

Te-^^Efe-^^'^W] dt = e-"^ (^ --)+-, (6.1) 

Jo Wss s) s 

r e-'' E \e-P(i-y'^m dt = e-"^ (^ --]+-. (6.2) 

Jo \Vs§ sj s 

In particular, for a = 

roo 1 

e-^*E[e-^^oW]dt = -p=. (6.3) 



Proof. It is sufficient to prove formula ()6.ip only; indeed, 

/ e-^*E[e-''(*-^"W)]dt= / e-^'*E[e^^"W]dt, (6.4) 

Jo Jo 

hence the left-hand side of (j6.2p can be computed using (16. ip by changing s to s and /3 to — /?, which 

amounts to interchanging the symbols s and s in (|6.ip . thus leading to formula (|6.2p . (Note that 

the right-hand side of (j6.4|) is well defined since Ya{t) ^ t and so e"'^**+'^)*E[e^^''(*)] ^ e"***.) 

Now, if a = then lo(i) = ^^i where Yq has the arcsine distribution with the density (II. Sp . 

hence the left-hand side of ()6.3p is reduced to 

/■oo /if* e~i^y \ 

/ e"^* - / dy dt. (6.5) 

Jo \^Jo VW^) ) 

The internal integral here can be interpreted as the convolution (/i * /2)(i) of the functions /i(i) = 
g-^t^-i/2 ^,^^ f2{t) = t~^'^, hence the Laplace transform (j6.5p reduces to the product 

1 /""e-^*e-^*t-i/2dt Te-^^t-i/^dt: ^^^^ ^^^^ ^ 



TT Jo Jo TTy/s + /3 y/s y/^S ' 

and the required formula (]6.3p follows. 

If a > then, noting that Qa/^iu) = tqa{ut) and using ()2.17p . ()2.18p and ()2.19p . we have 



E[,-/3y.W] = / q^i^u) du + / e-^y ( f'^^p^ ( -^ ) du ) dy. (6.6) 



t±- CIO 1 ■ 

— U \t — U 
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Interchanging the order of integration and making the substitution y = z{t — u), we can rewrite the 
second (iterated integral) term on the right-hand side of (j6.6p as 

J qa{u)(j e-^^'-^'>'p^{z)dz\du, 

which can be viewed as the convolution [qa * 'pp){t), where 

pp{t) := [ e-^''p^iz)dz = E[e-^^°W]. (6.7) 

Jo 

Returning to (j6.6p and applying the Laplace transform (with respect to the variable t), by the 
convolution theorem the left-hand side of (j6.ip can be expressed as 

■oo 

(6.8) 



/"OO / /"CXD \ /'CO /"OO 

/ e-^M / qa{u)du]dt+ e-'^qa{t)dtx e-'*pp{t)dt. 

Jo \Jt / Jo Jo 

Recall that, according to (I2.16p . 

poo 

/ e-^*g,(t)dt = e-'^v^, (6.9) 

Jo 

whence, integrating by parts and using (|2.16p . we obtain 



OO / roo \ 11 

e~"* / qaiu)du]dt= e^""^'. (6.10) 

\Jt J s s 



Furthermore, from (j6.7p and (|6.3p we have 



OO 1 

St - - ' - J- 



e"'*p;3(t)di = ^=. (6.11) 

Vss 

As a result, substituting expressions (|6.9p . (j6.10p and ()6.1ip into (j6.8p . we obtain formula (j6.ip . D 

Proof of Theorem 12.31 As T — )■ oo, we have ^ := {cT)~^x = (AT)^^'^(a-l-o(l)), whereas from (14. 7p 
it follows that k{s) ~ (2ATs)^/^ k{s) ~ (2Ars)i/2. Hence, from ()i30]) we obtain, for ^ ^ 0, a ^ 0, 

lim w${^, s) = -e»v^ ^ ^!"' + i = e-^v^ f^ --!+-> 

and similarly, for ^ ^ 0, a ^ 0, 

lim wUt s) = e-*^^ ^'"' ^ + - = e--^^ f^ - - V - • 



Comparing these results with Lemma 16. H by the continuity theorem for Laplace transforms 
|Fe71j we conclude that, for each t > 0, the distribution of the random variable trjj.^{x) (see (|2.ip 
and p.2p ) converges weakly, as T — )• oo, to the arcsine distribution (|1.3p if a = and to the 
distribution of either Y_a{t) if a < (see (j6.ip ) or t — Ya(t) if a > (see (j6.2p ). Speciahzed to the 
case t = 1, this readily gives the result of Theorem 12.31 Thus the proof is completed. D 

7. Proof of Theorem 12.41 

Similarly to SectionUl let us set wt{^, s) := (wj^(^, s), w^{^, s))' and rewrite equations (j3.1ip . ()3.12p 
in the matrix form (cf. (|4.ip . (|4.2p ) 



d^ 



AT{C,s;f)wT{^,s;f)-l, CeR\Df, (7.1) 



lim WT{C,s;f) = s-^l, lim WTi^,s; f) = s^H, (7.2) 

5— >— oo 5"^+°*^ 
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with the matrix (cf. ()4.3p ) 

Arii, s; f) := ATJi + {s + f3f{cn)) J2 , 

where Ji and J2 are defined in (I4.5p . Let us set 

6t{^, s) := wt{^, s; /) - wt{^, s; H), (7.3) 

where H is the Heaviside step function (cf. (j4.3p ). Owing to the properties of the solution wt{S,, s; f) 
(see the end of Section [3|), the function St{^,s) is bounded and continuous in ^ G M (for any fixed 
T, s > 0). Prom the relation (I7.3P and equations (I7.ip and (]7.2p . we obtain the differential equation 



dc 



AT{^,s;H)STi^,s) + foicT()wTiC,s), ^ eR\{DfU {0}), 



(7.4) 



where /q := f — H and Wt '■= /3J2Wt (for short), with the boundary conditions 

lim 5t(^,s) = 0. (7.5) 

5— !>±00 

More explicitly, equation (|7.4p splits into two equations on the negative and positive half- lines: 

dST{^,s) 






ATdTiC,s) + fo{cTOwT{C,s) 

At ST{ts) + fo{cm writs) 



e<o, 



(7.6) 
(7.7) 



where At = At{s) = XTJi + SJ2, At = At(s) = ATJi + s J2 (cf. g3D). 

By the variation of constants, equation (j7.6p is equivalent to the integral equation 



Sri^s) 



e^^^CT + 



.(i-y)AT , 



gv. -'■■-" fQ{cTy)wT{y,s)dy, ^ s^ 0, 



(7.8) 



where ct = cxis) = lim^^^o- ^t{£,,s) is a constant vector (for fixed T and s). By the exponential 
formula (j4.9p . equation (|7.8p takes the form 



Srits) 
where 



2 



(/ + K^UT)cT(s)+qJ(e,s 



2 



(/-K UT)cr(s)+qy(e,s) 



q$iC,s):={I±K-^AT) [ e^^yUcTy)wT{y,s)dy, ^ ^ 0. 

JO 



(7.9) 



(7.10) 



For fixed s and T, we have q^{t s) = e '^^ o{l) as ,^ — ?• —00. Indeed, via the change of variables 
z = y — ^ and applying Lebesgue's dominated convergence theorem, we see that 

0(1) / e-''^|/o(cT(z + 0)|dz = o(l), ^ ^ -00, 



-K(y-5) 



fo{cTy)wT{y,s)dy 



since wt and /o are bounded whereas /o(cT(z + ^)) — )■ for each z, according to the hypothesis of 
Theorem 12.41 Hence, due to the boundary condition (j7.5p at ^ = —00, equation ()7.9p implies 

e-€|(/_K-iAT)cT(s) + q^(^,s)}=o(l), ^ ^ -00. (7.11) 

Note that the expression in the curly brackets in (j7.1ip has a finite limit as S, — ?• —00, which then 
must vanish in order to extinguish the multiplier e~'^^ — t- 00, that is, 

{I - K^^AT)cT{s) = -q^{-oo,s). (7.12) 

Conversely, condition (j7.12p implies the limit (|7.1ip . since, by the I'Hopital rule, we have 

qTits)-q:p{-oo,s) _^ fo{cTC)wTits) 

~ (i — K At) = o(l), ^ — )• —00. 



aftg 



K 
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Analogous considerations applied to (|7.7p lead to the integral equation 

/■? 
ST{C,s)=e^^^CT+ e^^-y^^^McTy)wT{y,s)dy, ^ ^ 0, (7.13) 

Jo 

with ct = ct{s) = lim^_i.o+ ^riC, s), which, similarly to (j7.12p . implies the condition 

(/ + K"^iT)cT(s) = -gJ(+oo,s), (7.14) 

where k = k{s) and 

f< 
q^{ts):={I±k-'AT) e'^^yfo{cTy)wT{y,s)dy, ^ ^ 0. (7.15) 

Jo 

Moreover, since the function St{-,s) is continuous at ^ = 0, from formulas (j7.8p and (j7.13p we see 
that ct{s) = ct{s). Using this and subtracting (|7.12p from (|7.14p . we obtain 

ct{s) = [k'-^At + k~^At) q^{—oo,s) — qj,{+oo,s) . (7.16) 

Evaluating the matrix inverse in (j7.16p is facilitated by introducing the matrices (suggested by 
formulas (j4.6p ) 

and observing that 

' s/k\ 



This gives 



R-^AtK = K , 

\ k/s 



11 1 ~ 1 I ss \ 

K-^{K-^AT + K-^AT)K = {sk + sK)n^^{s), 7^T(s) := , (7.18) 



and, returning to ()7.16p . we finally get 

ct{s) = {sk + SK)~^KnT{s)K'^ q^{-oo,s) - q:^{+oo,s) . (7.19) 

In view of Theorem 12.31 and according to (17. Sp , to complete the proof of Theorem 12.41 we have 
to check that if ^yXT — T-aGMasT— )-oo then 5t{^, s) — )• 0. To this end suppose, for instance, 
that i ^Q and a ^ (the mirror case ^ ^ 0, a ^ is considered similarly). Note that, as T — )• oo, 

K ~ V2sAT, Ki = V2'sa + o{l), k;~Ut = ATk^^ Ji + 0(k~^). (7.20) 

Recall that the vectors qrp{£,,s), q^{S,,s) are defined in ()7.10p . (|7.15p . respectively. 
Lemma 7.1. For each s > 0, q^{—oo, s) = o(l) and g^ (+oo, s) = o(l) as T — )• cxd. 

Proof. Both q^ and q^ are considered similarly. For instance, using (|7.20p and making the change 
of variable z = Ky, we have 



fO 

\q^{-oo,s)\=0{l) e'\fo{cTK-h)\dz = o{l), T 

J —oo 



cx), (7.21) 



since, by the assumption of Theorem 12. 4[ cTk ^ ~ (2s) ^/^y^c^T/A — )■ oo, hence fo{cTHi ^z) — > 
for each 2 < 0, and we can apply Lebesgue's dominated convergence theorem. D 

Lemma 7.2. As T -)■ 00, if ot ■= S,VxT -^ a eR then, for each s > 0, q^{£,, s) — )• 0. 
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Proof. By the substitution y = ^z and with the help of asymptotic relations ()7.20p . we have 

Jo 

= 0{l)aT f \fo{zaTVc^T/X)\dz. (7.22) 

Jo 

Now, if ot — > a = then the right-hand side of (j7.22p vanishes in the limit as T — t- oo, since the 
function /q is bounded. If ay — ?• o 7^ then, similarly to the proof of Lemma 17. H the integral in 
()7.22p tends to zero thanks to Lebesgue's dominated convergence theorem. D 

Let us now return to equation (j7.9p . Using the identity ()7.12p and regrouping, we have 

St{^, s) = e'^^K^^ArCTis) - cosh«) ^^(-OO' «) + 2 ^"^ ^t(^' ^) + 2 ^"''^ ^t(^' ^) 

= 0(l)K-UTCr(s) + o(l), T^oo, (7.23) 

according to the second asymptotic relation in (|7.20p and Lemmas 17. II and 17.21 Further, substituting 
the expression ()7.19p for ct and using Lemma 17.11 and the last relation in (|7.20p , we obtain 

k-^AtCt = —, — ^ -(XTJi + 0(l))KnTK^^o(l), T ^ 00. (7.24) 

K{SK + SK) 

Li turn, using the expressions (I7.17P for the matrices K and K^^ and recalling the definition of the 
matrix TZt given in (j7.18p . it is easy to calculate 

KTZtK-^ = '^Ji + 0{l), T^oo. (7.25) 

Finally, combining (|7.24p and (|7.25p and noting that Jf = (see (03])), we have k~^AtCt = o(l) 
and hence, from (j7.23p . St{£,, s) = o(l) as required. This completes the proof of Theorem 12.41 

8. Concluding remarks 

We performed computer simulations to illustrate numerically the convergence to the arcsine law, 
as stated by Theorems 12.31 and 12.41 for the occupation time functional f?^(0;/) = T~^ L f{X^ )dt 
with various probing functions /. The simulation algorithm is easily implemented by virtue of the 
obvious decomposition 

f-T "■"! rn+i fT~<Jn 

T4io,f)^ fixt)dt = y2 f{x^,+i-iyct)dt+ f{x^„+{-irct)dt, 

Jo j^Q Jo Jo 

where (tj) is a sequence of independent random times with exponential distribution each (with 
parameter A), and ctj := ti + • • • + r, are the successive reversal times of the telegraph motion; the 
threshold value n is determined by the condition (t„ ^ T < (T„ + t„+i. 

Throughout the simulations, we used the standardized parameters c = 1, A = 1, and plotted 
histograms of the sample values of r]j, (0; /) based on A^ = 10, 000 runs of the telegraph process. To be 
specific, we simulated the plus-version of the process, X^ (i.e., with positive initial velocity), leading 
to histograms slightly skewed to the right, especially at moderate times T. No formal goodness-of-fit 
tests were applied, but the histograms in Figures [U and [2] below clearly demonstrate the developing 
[/-shape characteristic of the arcsine distribution, however with the speed of such a convergence 
apparently depending on the function / involved (and, of course, on the observation time T used). 

We start with the "canonical" case where the Heaviside function H{x) = l(o,oo)(2;) plays the 
role of the probing function /. Simulated values of 7]^{0;H) were obtained over the observation 
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time T = 1 000. The histogram plotted in Figure [1^ shows a very good fit of the data to the 
theoretical arcsine density (rescaled according with the chosen representation of the histogram). 
As already mentioned, the noticeable difference between the highest columns at the left and right 
edges may be attributed to asymmetry of the process X^ . More precisely, the proportion of the 
sample values of r]^{0;H) falling, say, in the first box, Ai (from to 0.01) and the last box, Aioo 
(from 0.99 to 1) is given by 510 and 750, respectively, yielding the relative frequencies 510/10,000 = 
0.051 and 750/10,000 = 0.075. The corresponding limiting probabilities, computed from the arcsine 
distribution ()1.3|) . equal 0.064 for both Ai and Aioo (here and below, we give numerical values to 
two significant figures). This discrepancy can be quantified using the exact theoretical distribution 
of r/^(0; H) obtained in Theorem 12. II (see formula (|2.8p with T = 1000), giving the probability 0.052 
for Ai and 0.077 for Aioo, where the latter includes the atom 2(/9t(1) = 0.025. For comparison, with 
a tenfold observation time T = 10000, these probabilities become 0.060 and 0.068, respectively, with 
the atom much reduced, 0.008. It is also worth mentioning that, as indicated by these results, the fit 
with the limiting arcsine distribution would be much better for the "symmetric" version r}T{0;H) 
corresponding to the telegraph process Xt (see ()1.5p ). 
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Figure 1: Histograms for the occupation time functional 7?J(0;/) with (a) the Heaviside step 
function f = H and (b) the function f{x) = 7r~^arctan3; + i. The parameters of the telegraph 
process X^ are standardized to c = 1 and A = 1. Both histograms are obtained with N = 10,000 
simulations, each over the observation time T = 1000. The length of each box on the histogram is 
A = 0.01. The red solid curve represents the scaled arcsine density (i.e., multiplied by A^A = 100). 



The long-term prediction contained in a more general Theorem 12.41 was verified by computer 
simulations for the functional 'i]^{0]f) with the probing function f{x) = 7r^^arctanx+ ^. The new 
histogram plot (see Figure [Tja), obtained with the same values of c. A, T and N, is qualitatively 
similar to that on Figured^, including a small right bias, but convergence to the arcsine distribution 
becomes slower, apparently due to additional time needed for the process to explore the limiting 
values f± of the function / at ±00, which eventually determine the distributional limit. 

Incidentally, this observation helps to understand the difference between the sets of hypotheses in 
Theorems 12.31 and 12.41 : indeed, the additional condition of Theorem 12. 4^ requiring that (?T/X — )• 00 
as T — )• 00, guarantees a sufficient mobility of the telegraph process needed to gauge the limits f± 
available only at remote distances from the origin. In contrast, if the function / is reduced to the 
Heaviside step function H^ the limiting values H^ = 0, if+ = 1 are encountered by the process 
straight away, so no extra mobility is needed. 

Let us point out that the asymptotic conditions imposed in Theorem 12.41 on the function / 
involved in the occupation functional 7]j,{x; f) are rather strong, assuming the existence of the limits 



19 



Leonid Bogachev and Nikita Ratanov 



350 
300 



100 
50 




0.4 0.6 0.8 1 




Figure 2: Histograms for the functional ?7^(0; /) with the probing function f{x) = 7r~^arctanx + 
cos X + 2 ■ The parameters of the telegraph process are as in Figure [H with the same number of 
runs A^ = 10,000 and the observation time (a) T = 1000 or (b) T = 10000. Compare with Figure 
[J and note the improved quality of fit to the hypothetical arcsine distribution (red curve) on the 
right plot as compared to the left one. 



lima;_5.-|-oo fix) = f± (with /_ ^ /+) • This is in contrast with the paper by Khasminskii |Kh99j 
mentioned in the Introduction, where the function / is only assumed to be Cesaro (C, l)-summable, 
i.e., subject to a weaker condition lim2,_s.-|-oo ^~^ Jq /("") du = f± (cf. ()1.4p ). Unfortunately, we were 
unable to reach the same level of generality. In particular, our proofs of formulas (|7.12p . (j7.14p and 
the key Lemmas 17.11 and 17.21 (see Section [7]) are heavily based on the existence of the limits f± . 

However, we conjecture that Theorem 12.41 does hold under the weaker condition of Cesaro (C, 1)- 
summability of the probing function /. To verify this claim numerically, we carried out computer 
simulations for the distribution of ?7j(0; /) with f{x) = vr^^arctanx + cosx + ^ ■ Figure [2]a shows 
the simulated histogram with the old values T = 1000 and A^ = 10,000, which reveals a bimodal 
distribution but not quite well fit to the hypothetical arcsine limit; in particular, there are noticeable 
"parasite" shoulders outside the interval [0,1], which are indeed possible because the function / 
may take values less than and bigger than 1. However, the fit with the arcsine shape significantly 
improves under longer observations, T = 10000 (see Figure [2]b). In particular, the high modes at 
the edges are better pronounced, while the shoulders outside [0, 1] are considerably reduced. 



Appendix A. Probabilistic proof of Theorem 12.21 

Let us recall some information related to the first-passage problem for the telegraph process X^ . 
For X < 0, let T_^. := min{t ^ : X^ = —x} (with the convention that inf := +oo) be the hitting 
time of point —x > by the process X^ (starting from the origin, X^ = 0). If we set Tq := {—x)/c, 
then the distribution of 1-_^ is concentrated on [To,oo) and is given by (see |Pi9H §0.5, pp. 12-13] 
and also |Or951 Theorem 4.1, p. 18]) 

P{Tl^ G dt} = e-^^«(5To(dt) + Ql^(t) dt, ¥{1Z^ e dt} = QZ^{t) dt, (A.l) 

where the densities Q_^ are defined exactly by equations (j2.13p . (j2.14p . 

Consider the two-dimensional Markov process [X^ ,V^ ) , where Xf. is the (conditional) telegraph 
process (jl.6p (i.e., with the initial velocity Vq = itc, respectively), and V^ = dXj /dt = ibc(— 1)^* 
is the corresponding velocity process driven by an dependent Poisson process Nt which determines 
the reverse instants of the motion X^ (see (jl.6p ). It is obvious that T.^. is a stopping time for the 
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process (X. , K ). Also note that K I ~± = +c (a.s), since the first passage through point —x > 

by the process X^ , starting from the origin, with probabihty 1 can only occur from left to right, 
that is, with positive velocity. Hence, conditioning on the hitting time of the origin starting from 
X < (which, of course, has the same distribution as T^^,) and using the strong Markov property 
of the joint process {X^ , VJ ), we have, for each y G [0, 1 — Tq/T], 

P{r/±(rc) G dy] = ¥{%% > T}6o{dy) + E[P{r/±(a;) G dy, To^lt,^T\ T^J] 

/ P{T±, G du] 5o{dy) + / P{T±, G du] P{(1 - u/T) r/+„ JO) G dy}. (A.2) 

\Jt J jTq 

Here, the first integral represents the case where the telegraph process X^ does not reach the origin 

before time T and, therefore, never enters the positive half-line (thus contributing to the atom 

(^o(dy)), while the second integral (where integration is taken with respect to du) accounts for the 

first passage event (at time instant u G [Tq, (1 — y)T]), so that the telegraph process, restarted from 

the origin (with the initial velocity +c), has to spend on the positive half- line the required time 

Tdy during the remaining travel time T — u. 

In view of (|A.ip together with ()2.13p and p.l4p . and due to equation p.Sp which provides the 

distribution of ?7j_^(0), formula ()A.2p furnishes an explicit representation of the distribution of 

T]rp (x) . More explicitly, on account of the atom in (jA.ip , the right-hand side of (|A.2p specializes to 

(1-?^)^ , r . dy 



Q^,iu)duJ6o{dy) + 4idy) + J^ Q±J^) P |^+_J0) G ^— ^ j dn, (A.S) 

where //^(dy) := and 

/i+(dy) := e-^^« P Jt^+.^JO) G J^rj^] ■ (A.4) 

Using (j2.8p . for any u G [Tq, (1 — y)r] we have 

Substituting (|A.5p (with u = Tq) into (jA.4p readily gives (|2.1ip . while the last term on the right-hand 
side of (TOl) is reduced to (cf. (l230l) ) 

2rQ± j(i - y)T) My) ^y + £ ' Q-M {^T-u (y3^) Y~^) *^''' 

where the contribution of the atom (5i_u/j'(dy) from (jA.SP is easily computed via the obvious 
symbolic formula (Ji-n/T(dy) du = T(5(i„y)7-(d«) dy. Indeed, for any test functions F{y) and G{u) 
we have, by changing the order of integration, 

l-To/T r(l~y)T I'T /-l-u/T 

F{y) Giu)6,_^/T{dy)du= G{u)du F{y) d^.^/ridy) 

Jto Jto Jo 

[■T rl^To/T / i-{l~y)T \ 

= G{u)F{l-u/T)du = T F{y)[ G{u) 6(^,^y)T{du) ] dy. 

Jtq Jo \Jto J 

Appendix B. Probabilistic proof of Theorem 12.31 

Making the substitution t = Tu and using that H{ax) = H{x) for any a > 0, we can rewrite 
formula (j2.ip as 

-1 ri 



7?^(x)= I H{x + X^Jdu= I H{3: + Xt)du, (B.l) 

Jo Jo 
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where x := {c^T /\)~'^''^x, and X^ := {c^T/X)~^''^X^^ [u ^ 0) is another telegraph process with 
rescaled parameters A := AT — )• oo, c := (XT)^'"^ — )• oo (T — )• oo). By Theorem II .1^ the process 
{X^, ^ u ^ 1) converges weakly to a standard Brownian motion {B^, ^ u ^ 1). Hence, if 
X — )• a as T — )■ oo (cf. the hypotheses of Theorem I2.3P then from (jB.ip we immediately obtain the 
convergence in distribution, as T — )• oo, 



Jo 



According to (jl.ip and (|1.2p . the random variable f)i(O) has the arcsine distribution, which proves 
Theorem 12.31 for a = 0. For a < (so that — a > 0), let T_a := min{t ^ : Bt = —a} be the hitting 
time of the point —a by the Brownian motion Bt starting from the origin {Bq = 0). As is well known 
since P. Levy's paper [LiiOl Theoreme 2, p. 294] (see also ^M7^ §1.7, p. 26] or [E^TTI §VL2(e), 
pp. 174-175]), the random variable r_a has probability density q-a{-) defined in (j2.15p . Note that 
T-a is a stopping time (with respect to the natural filtration J^t '■= (^{Bg, ^ s ^ t}). Conditioning 
on T_a (when a + Bt-_^ = 0) and using the strong Markov property, we obtain, for any y G [0, 1], 

P{(li(a)Gdy}=P{r_, >l}<5o(d?/)+ / %_a(n) P{(1 - n) f)i„„(0) G dy} du 

Jo 

q-a{u) du ] 6oidy) + ( / -p pas ( :; I dn ) dy, 



— u \1 — U^ 

which coincides with (j2.17p (for Y_a) in view of (|2.18p and (|2.19p . Finally, the case a > easily 
follows by noting the obvious symmetry relation f)i(a) = 1 — fii(— a) (cf. (12. Sp ). 
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